Elliptic flow of A hyperons in Pb+Pb collisions at 158A GeV 
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The elliptic flow of A hyperons has been measured by the NA49 collaboration at the CERN-SPS in 
semi-central Pb+Pb collisions at 158^4 GeV. The standard method of correlating particles with the 
event plane was used. Measurements of vi near mid-rapidity are reported as a function of rapidity, 
centrality and transverse momentum. Elliptic flow of A particles increases both with the impact 
parameter and with the transverse momentum. It is compared with V2 for pions and protons as 
well as with model calculations. The observation of significant elliptic flow and its mass dependence 
suggest strong collective behavior of the matter produced in collisions of heavy nuclei already at the 
SPS. Scaling properties of elliptic flow of different particle species have been tested at 158^4 GeV. 
The limited pt range of the data does not allow for a decisive test of the coalescence model. 
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I. INTRODUCTION 

Elliptic flow in relativistic nuclear interactions has its 
origin in the spatial anisotropy of the initial reaction 
volume in non-central collisions and in particle rescat- 
terings in the evolving system which convert the spatial 
anisotropy into a momentum anisotropy [l| . The spatial 
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anisotropy decreases rapidly because of the fast expan- 
sion of the system Q making the momentum anisotropy 
measured at the end of this evolution strongly depen- 
dent on the matter properties and the equation of state 
(EoS) at the early stage [1, 0|- Comparison of measured 
anisotropics with hydrodynamic model calculations pro- 
vide an important test of the degree of thermalisation in 
the produced particle system at the early stage. Flow of 
heavy particles is affected more strongly by changes in 
the EoS than flow of pions [4|, [5j, [6J . 

The anisotropic flow parameters measured to date at 
SPS and lower energies are mainly those of pions and pro- 
tons [7J,|8j]. RHIC experiments have measured elliptic flow 
for many particle species [9| , including hyperons @, [l(| • 
In these experiments, the rapid rise of elliptic flow with 
transverse momentum (px) up to 1.5 GeV/c and its 
particle-mass dependence are well reproduced by hydro- 
dynamic models [9J . At higher values of pr quark number 
scaling of elliptic flow has been observed at RHIC [13, [Hj] 
which indicates that the quark coalescence mechanism 
dominates hadron production in the intermediate px re- 
gion [12j| . In order to test the validity of this scenario 
at SPS energies we have extended elliptic flow measure- 
ments in 1584 GeV Pb+Pb (y/s^ = 17.2 GeV) colli- 
sions to A hyperons. In this paper, results on A elliptic 
flow as a function of the center of mass rapidity (y) and 
transverse momentum will be presented and compared to 
model calculations. 



II. ANALYSIS 

The main components of the NA49 detector [l3j] are 
four large-volume Time Projection Chambers (TPCs) 
for tracking and particle identification by energy loss 
(dE/dx) measurement with a resolution of 3—6%. The 
TPC system consists of two vertex chambers inside the 
spectrometer magnets and two main chambers placed be- 
hind the magnets at both sides of the beam. Downstream 
of the TPCs a veto calorimeter detects projectile spec- 
tators and is used for triggering and centrality selection. 
The data sample consists of 3 x 10 6 semi-central Pb+Pb 
events after online trigger selection of the 23.5% most 
central collisions. The events were divided into three dif- 
ferent centrality bins, which correspond to the first three 
bins used in a previous analysis (see Table 1 in Q) and 
are defined by centrality ranges 0-5% (bin 1), 5-12.5% 
(bin 2), and 12.5-23.5% (bin 3). Many model predic- 
tions are published for impact parameter ranges similar 
to those of our centrality classes, which are: 0-3.4 fm (bin 
1), 3.4-5.3 fm (bin 2), and 5.3-7.4 fm (bin 3) for our cen- 
trality classes. The measurement in the centrality range 
c/ctot = 5—23.5% (called mid-central) is obtained by 
averaging the results of bins 2 plus 3 with weights cor- 
responding to the fractions of the total cross section in 
these bins. 

The A hyperon candidates were selected from the sam- 
ple of V°-track configurations consisting of oppositely 



charged particles, which include the A decays into pro- 
ton and 7T _ (branching ratio 63.9%). The identification 
method [l4[ relies on the evaluation of the invariant mass 
distribution and is enhanced by daughter particle identi- 
fication applying a cut in dE/dx around the expectation 
value derived from a Bethe-Bloch parametrization. The 
extracted A candidates have a background contamination 
of 5-9% in the p-7r invariant mass window 1.108—1.124 
GeV/c 2 , depending on centrality. The yields of A hyper- 
ons are obtained by counting the number of entries in the 
invariant mass peak above the estimated background as 
a function of the azimuthal angle of A candidates with 
respect to the event plane angle $2EP- The background 
is estimated from a fit of invariant mass spectra to the 
sum of a Lorentz distribution and a polynomial. It is 
subtracted in 16 bins of tf> — < I > 2EP- The acceptance of 
A hyperons covers the range 0.4 < px ^ 4 GeV/c and 
— 1-5 J$ V ^$ 1.0 but the detection efficiency strongly de- 
pends on azimuthal angle, px and y. Thus we have to 
introduce differential corrections to avoid biases when av- 
eraging V2 over rapidity and transverse momentum. Mul- 
tiplicative factors were introduced for every A particle to 
correct the A yields for detector and reconstruction effi- 
ciency. They were determined as ratios of published A 
yields jla |. parametrized by the Blast Wave model, to 
the measured raw A yields. 

The elliptic flow analysis is based on the standard pro- 
cedure outlined in [1, [l6| to reconstruct the event plane 
for each event and the corrections for the event plane 
resolution. The event plane is an experimental estimator 
of the true reaction plane and is calculated from the az- 
imuthal distribution of primary charged it mesons. Iden- 
tification of pions is based on dE / dx measurements in the 
TPCs. To avoid possible auto-correlations, tracks asso- 
ciated with A candidates are excluded from the event 
plane calculation. The method to determine the event 
plane azimuthal angle $2EP uses the elliptic flow of pi- 
ons, according to the formula: 

N 

X 2 = £> T [cos(2<^) - <cos(20)>], 

z=l 
N 

Y2 =X>TMn(2^)-<sm(20))], (1) 

z=l 

* 2 EP 4 tarl (l)' 

where X 2 , Y% are the components of the event plane flow 
vector Q2 and the sums run over accepted charged pion 
tracks in an event. The acceptance correction is based 
on the recentering method Q which consists of subtract- 
ing in Eq. 1 the mean values (cos(2</>)) and (sin(20)). 
These mean values are calculated in bins of px and ra- 
pidity for all pions in those events which contain at least 
one A hyperon candidate. The values were stored in a 
3-dimensional matrix of 20 px intervals, 50 rapidity in- 
tervals, and eight centrality bins. A second level accep- 
tance correction is done by using mixed events. We used 
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FIG. 1: (color online). Azimuthal distributions of A hyperons 
with respect to the event plane for real events (solid symbols) 
and mixed-events (open symbols) in two centrality bins. The 
curves are Fourier expansion fits (see text Eq. (2)). 



10 mixed events for each real event. Particles for mixed 
events are randomly selected from different events in the 
same centrality bin with at least one A hyperon. The final 
angular distributions are obtained by dividing the real A 
angular distribution by the mixed event distribution to 
remove the acceptance correlations remaining after re- 
centering. The corrected A azimuthal distributions are 
then fitted with a truncated Fourier series: 

dN 



d{(j> - $2EP) 

const x (f 



+ uf s cos[2(^-$ 2EP )] 
+ < bs cos[4(0-$ 2EP )]). 



(2) 



The elliptic flow v-i is evaluated by dividing the observed 



anisotropy uf hs by the event plane resolution R: 



V2 



obs 
H 

R ' 



(3) 



The resolution, 



R = J2(cos[2($° 



2EP 



*Lp)]>. 



is calculated from the correlation of two planes 
(^EP' ^ep) f° r random sub-events with equal multiplic- 
ity. The results are R = 0.27, 0.34 and 0.40 for centrality 



bins 1, 2 and 3, respectively. The total errors in Figs. 2- 
4 are given by the quadratic sum of contributions from 
the statistical error of the signal, the uncertainty due to 
the background subtraction, the mixed event correction 
and event plane resolution. The observed hexadecupole 
anisotropy v\ hs is consistent with zero within statistical 
errors for all centrality bins. 



III. RESULTS 

The final statistics in our sample consists of about 10 6 
A candidates. This allows flow analysis for several rapid- 
ity and pt bins. Two sample azimuthal distributions of A 
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FIG. 2: (color online) . Elliptic flow of A hyperons as a function 
of rapidity (top) and px (bottom). The open points in the 
top graph have been reflected about midrapidity. 



hyperons with respect to the estimated reaction plane for 
real and mixed events are shown in Fig. f. The curves 
represent results of fits with the truncated Fourier se- 
ries Eq. (2). The distributions exhibit a strong corre- 
lation for real events (full symbols and curves). As ex- 
pected, no correlation is observed for mixed-events (open 
symbols, dashed curves). The pt averaged elliptic flow 
is obtained from all identified A hyperons without pt 
cuts. It exhibits no significant dependence on rapidity as 
shown in Fig. 2(top). The absence of a rapidity depen- 
dence of i>2(y) was also observed for protons (see Fig. 6 
in Ref. [H) in mid-central events. We use the A sample 
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FIG. 3: (color online). Elliptic flow of A hyperons as a func- 
tion of pt from mid-central (top) and central (bottom) events 
measured by the STAR (open symbols) and NA49 (solid sym- 
bols) experiments. Curves are hydrodynamical model predic- 
tions at RHIC energy for the two different centrality bins. 



from the full rapidity range of the data in Fig. 2(top) 
for the study of vi as a function of px. The px de- 
pendence of rapidity-averaged A elliptic flow is shown 
in the bottom plot of Fig. 2 for two centrality ranges. 
The V2 parameter significantly increases with transverse 
momentum, the rise being stronger for more peripheral 
events. Fig. 3 shows a comparison of V2(px) of A hyper- 
ons for mid-central and central (0-5%) events measured 
by the NA49 and STAR experiments [l7j]. The NA45 
collaboration recently also presented preliminary results 
in Pb+Au collisions at the top SPS energy [18| which 
agree well with the NA49 results (not shown). For mid- 
central collisions at SPS energy the elliptic flow grows 
linearly with px up to ^2 GeV/c, but the increase is 
steeper at RHIC than at SPS energy. It should be noted 
that RHIC mid-central data have been measured in the 
centrality range ct/oxot = 5—30% while SPS events are 
somewhat more central. The effect of different central- 
ity ranges has been estimated by hydrodynamic calcula- 
tions H3| at RHIC energy for the slightly different 
centrality bins of NA49 and STAR. As shown by the cor- 
responding curves in Fig. 3 this explains only partially 
the difference between both measurements. For central 
collisions NA49 and STAR results do not differ signifi- 
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FIG. 4: (color online). Elliptic flow for charged pions (cir- 
cles) and protons (triangles) as a function of px from mid- 
central events measured by the STAR [t| (open symbols) and 
NA49 8] (solid symbols) experiments. 



cantly (Fig. 3 bottom). A comparison of vi of pions and 
protons obtained by NA49 [|] and STAR for mid- 
central collisions is shown in Fig. 4. The NA49 values 
were obtained as the cross section weighted averages of 
the measurements published in Sj for the appropriate 
centrality range. As observed for A hyperons, v<x of pions 
also rises faster with px at RHIC than at SPS. Protons 
seem to exhibit a similar trend, but the limited px range 
and the larger measurement errors do not allow a firm 
conclusion. 

A comparison of t>2(px) for pions, protons and A hy- 
perons as measured by the NA49 experiment in mid- 
central events is displayed in Fig. 5. The elliptic flow 
grows linearly with px for all particle species but the rise 
for pions starts from px close to zero while for protons 
and A particles it starts from px ~ 0.5 GeV/c. The el- 
liptic flow for pions is significantly larger than that for 
heavier particles although at px ~ 2 GeV/c the flow 
becomes similar for all particle species. Data are re- 
produced by blast wave fits @, EH (dashed curves m 
Fig. 5) with the following parameters: freeze-out tem- 
perature T — 92 MeV, mean transverse expansion rapid- 
ity of the shell po = 0.82, its second harmonic azimuthal 
modulation amplitude p a — 0.021, and a spatial eccen- 
tricity parameter S2 — 0.033. The values of freeze-out 
temperature and expansion rapidity are consistent with 
those obtained from mx spectra and Bose-Einstein cor- 
relations [22j . In Fig. 5 the measured values of V2 are also 
compared to hydrodynamical model calculations [23| as- 
suming a first-order phase transition to a QGP at the 
critical temperature T c = 165 MeV. The initial condi- 
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plain the quark number scaling observed at RHIC for 
t>2 in the intermediate px region. Fig. 6 shows the scal- 
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FIG. 5: (color online). Elliptic flow for charged pions (circles), 
protons (triangles) and A hyperons (squares) as a function of 
px from 158A GeV Pb+Pb mid-central events measured by 
the NA49 experiment. Curves are blast wave fits and hydro- 
dynamic model predictions for ^/^nn = 17-2 GeV. 



tions employed for the hydrodynamical calculations at 
SPS arc: initial energy density eo = 9.0 GeV/fm 3 , baryon 
density rif, = 1.1 fm~ 3 , thermalization time To = 0.8 fm/c. 
The initial conditions were fixed as in [2H by requiring 
a good fit to the px spectra of protons and negatively 
charged particles in central Pb+Pb collisions at the SPS. 
With the freeze-out temperature Tf = 120 MeV tuned 
to reproduce particle spectra, the model calculations sig- 
nificantly overestimate the SPS results for semi-central 
collisions (full curves in Fig. 5) in contrast to predic- 
tions at RHIC energy which agree quite well with data 
for pt < 2 GeV/c [17| (see Fig. 3). The discrepancy at 
SPS may indicate a lack of complete thermalisation or a 
viscosity effect. However, the model reproduces qualita- 
tively the characteristic hadron-mass ordering of elliptic 
flow. Thus the data support the hypothesis of early de- 
velopment of collectivity. On the other hand, the magni- 
tude and trends of elliptic flow, measured by the second 
Fourier coefficient t>2 is found to be underpredicted by a 
hadronic cascade model [25] . A more comprehensive de- 
scription of experimental results was found by coupling a 
hadronic rescattering phase to the hydrodynamical evo- 
lution and hadronisation This model can reproduce 
rax spectra and elliptic flow both at top SPS energy and 
RHIC consistently, although predictions of V2 for exactly 
the centrality range of the present analysis are not pub- 
lished in the literature. One should note, however, that 
none of the hydrodynamical models have yet been able 
to describe Bose-Einstein correlations successfully. 

Quark coalescence models [UJ have been used to ex- 
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FIG. 6: (color online). Elliptic flow for charged pions (circles), 
protons (triangles) and A hyperons (squares) scaled by the 
number of constituent quarks V2 /n as a function of px l n (top) 
and KE T /n (bottom) from 158A GeV Pb+Pb mid-central 
events measured by the NA49 experiment. 

ing behavior of the vi measurements of NA49 at the SPS. 
The values of v 2 shown in Fig. 5 were divided by the num- 
ber of constituent quarks n = 3 for baryons and n = 2 
for mesons. 

When plotting v^jn versus px /n approximate scaling 
was observed at RHIC @,[ll| except for pions. The NA49 
measurements (Fig. 6, top) are consistent with this result 
in the px range covered by the data. The deviation of 
the pions from the universal curve has been attributed 
to the Goldstone nature of the pion (its mass is smaller 
than the sum of masses of its constituent quarks) or to 
the effect of resonance decays [26j]. 
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The variable KE^ = mi — m, where ttzt is the trans- 
verse mass and m the rest mass of the particle, was 
proposed in [13] as an alternative to px since pressure 
gradients which give rise to azimuthal asymmetry may 
naturally lead to collective transverse energy of pro- 
duced particles. Good scaling for all particle species is 
seen when vijn is plotted versus KEt/u in the range 
KE T /n < 0.8 GeV covered by the SPS data (Fig. 6, bot- 
tom) . This behaviour was first observed at RHIC and in- 
terpreted as a result of hydrodynamic evolution [27] . The 
data of NA49 do not reach higher transverse momenta at 
which a decisive test of coalescence models would be pos- 
sible. 



IV. SUMMARY 

In summary, we report the first measurement of the 
anisotropic flow parameter V2 for A particles from Pb+Pb 
collisions at ^/snn = 17.2 GeV. Elliptic flow of A hy- 
perons exhibits no significant dependence on rapidity for 
— 1-5 ~ V ~ 1-0- It rises linearly with px and is smaller 
than V2 for pions. Both features are quantitatively re- 
produced by the Blast Wave parametrization but only 
qualitatively by the hydrodynamic model. The increase 
of v 2 with px is weaker at SPS than at RHIC energy. The 
observation of significant elliptic flow and its mass depen- 



dence suggest strong collective behaviour of the matter 
produced in collisions of heavy nuclei already at the SPS. 
Hydrodynamic models with a deconfinement phase tran- 
sition and a microscopic freezeout treatment appear to 
provide a consistent description of V2 and mx spectra 
at both top SPS and RHIC energies Quark number 
scaling of elliptic flow (vz/n versus KE^/n) was shown 
to hold also at the SPS. However, the limited px reach 
of the data does not allow a decisive test of the quark 
coalescence hypothesis. 
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